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摘要 : 明 安 图 射电 频谱 日 像 仪 (Mingantu Ultrawide SpEctral Radioheliograph，MUSER ) 是 
新 一 代 具 有 高 时 间 、 高 空间 、 高 频率 分 辨 率 的 太阳 专用 射电 望远镜 ， 采 用 综合 孔径 原理 成 
像 ， 所 以 幅度 和 相位 是 决定 最 后 成 图 质量 的 关键 因素 。 天 线 的 相位 方向 图 会 影响 日 像 仪 输出 
的 幅度 和 相位 ， 根 据 日 像 仪 的 馈 源 设计 和 综合 孔径 原理 ， 针 对 明 安 图 射电 频谱 日 像 仪 天 线 数 
目 多 ， 且 为 户外 环境 , 根据 天 文 观测 须 经 常 测试 天 线性 能 的 特点 ， 给 出 了 基于 相关 结果 测量 
日 像 仪 天 线 相 位 方向 图 的 方法 ， 该 方法 可 以 直接 通过 日 像 仪 的 相关 输出 结果 高 效 准确 地 得 到 
天 线 的 相位 方向 图 。 对 MUSER-I 天 线 的 相位 方向 图 进行 了 测量 和 分 析 ， 同 时 分 析 了 天 线 相 
位 方向 图 对 日 像 仪 成 像 的 影响 ， 为 得 到 高 质量 的 太阳 图 像 提 供 了 参考 和 保障 。 
关键 词 : 相位 方向 图 ; 相关 ; 明 安 图 射电 频谱 日 像 仪 
中 图 分 类 号 : P164 ”文献 标识 码 : A MBAS: 1672-7673(2016)03-0293-07 


天 线 的 方向 图 是 表征 天 线 辐射 特性 与 空间 角度 的 图 形 ， 如 幅度 、 极 化 、 相 位 等 。 通 常 是 以 天 线 的 
T= 相位 中 心 为 球 心 的 一 个 三 维 空间 图 形 ， 在 半径 为 的 足够 大 球面 上 ， 依 次 逐 点 测量 天 线 的 辐射 特性 给 
es 制 而 成 。 根 据 测 量 的 物理 量 不 同 ， 可 以 得 到 不 同 的 方向 图 。 如 测量 场 强 振幅 ， 得 到 场 强 方向 图 ; 测量 
O 功率 ， 得 到 功率 方向 图 ; WAAL, TARAN E 。 

天 线 的 相位 特性 包含 许多 非常 重要 的 信息 ， 特 别 是 随 着 相 控 阵 和 综合 孔径 技术 的 发 展 ， 天 线 相位 
测量 变 的 越 来 越 重要 。 天 线 相位 方向 图 测量 的 方法 同 幅度 方向 图 的 测量 方法 类 似 ， 所 不 同 的 是 要 测量 
待 测 天 线 在 远 区 球面 上 辐射 场 的 相位 。 通 常 测量 原理 与 一 般 的 相位 测量 类 似 ， 即 将 待 测 信 号 与 基准 信 
号 比较 得 出 相对 相位 值 ， 如 图 1°29), 
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图 1 传统 相位 方向 图 测试 框图 
Fig.1 Ablock diagram of measuring the phase pattern traditionally 
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但 是 明 安 图 射电 频谱 日 像 仪 天 线 阵 位 于 内 蒙古 大 草原 上 ， 为 野外 环境 ， 并 且 天 线 数目 众多 ” ， 
根据 天 文 观测 的 需求 须 经 常 测试 天 线 的 性 能 ， 希 望 能 够 快速 准确 地 得 到 天 线 的 相关 参数 。 因 此 本 文 给 
出 了 一 种 基于 相关 测量 天 线 阵 天 线 相位 方向 图 的 方法 。 首 先 ， 根 据 综合 孔径 原理 ， 推 导出 阵列 天 线 相 
位 与 可 见 度 函数 的 关系 ， 然 后 根据 明 安 图 射电 频谱 日 像 仪 的 相关 设计 ， 给 出 基于 相关 结果 测量 日 像 仪 
天 线 相位 方向 图 的 方法 。 并 对 MUSER-I 天 线 的 相位 方向 图 进行 了 测量 和 分 析 ， 同 时 分 析 了 天 线 相位 
方向 图 对 日 像 仪 成 像 的 影响 。 


1 明 安 图 财 电 频谱 日 像 仪 简 介 


明 安 图 射电 频谱 日 像 仪 是 一 台 对 太阳 进行 观测 的 专用 射电 望远镜 阵列 ， 在 厘米 -分 米 波段 同时 以 
高 时 间 、 高 空间 、 高 频率 分 辨 率 对 太阳 成 像 。 明 安 图 射电 频谱 日 像 仪 在 研究 瞬间 高 能 现象 、 日 园 磁 场 
和 太阳 大 气 结 构 、 太 阳 焰 斑 与 日 早 物 质 抛 射 的 源 区 特性 等 方面 发 挥 重要 作用 ， 从 而 使 人 类 更 好 地 了 人 解 
日 园 大 气 的 动力 学 过 程 "" 。 

明 安 图 射电 频谱 日 像 仪 建设 在 距离 北京 400 km 的 国家 天 文 台 明 安 图 观测 站 ， 根 据 综合 孔径 原理 ， 
由 低频 阵 和 高 频 阵 两 部 分 组 成 。 其 中 低频 阵 由 40 面 4.5 m 的 天 线 组 成 ,观测 频率 为 400 MHz 到 
2 GHz; 高 频 阵 由 60 面 2. 0 m 的 天 线 组 成 ， 观 测 频率 为 2 GHz 到 15 GHz。 所 有 的 100 面 天 线 分 布 在 三 
条 螺旋 辟 上 ， 高 低频 天 线 将 接收 的 太阳 射电 信号 通过 等 长 的 光纤 传送 到 室内 模拟 接收 机 ， 通 过 数字 接 
收 机 的 数字 化 处 理 到 达 复 相关 融 件 ， 对 所 有 通道 进行 复 相 关 计 算得 到 复 可 见 度 函数 ， 最 后 经 过 成 像 处 
理 得 到 太阳 图 像 '“ 。 

现 阶 段 明 安 图 射电 频谱 日 像 仪 已 经 完成 高 低频 阵 的 建设 工作 ， 系 统 正 处 于 调试 和 试 观测 过 程 中 ， 
最 新 的 主要 性 能 指标 如 表 1 7 。 

表 1 明 安 图 射电 频谱 日 像 仪 主 要 系统 指标 
Table 1 MUSER Specifications 


MUSER-I MUSER-II 
Frequency range 400 MHz-2 GHz 2 GHz-15 GHz 
Number of antennas 40 60 
Number of baselines 780 1770 
Antenna size 4.5m 2.0m 
Baseline range ~8 m-~3 000 m =4m-~3 000 m 
Polarization Dual circular polarization Dual circular polarization 
Dynamic range =25 dB =25 dB 


2 ”相位 方向 图 测量 


2.1 测量 原理 

明 安 图 射电 频谱 日 像 仪 是 采用 综合 孔径 技术 成 像 的 干涉 仪 系 统 ， 射 电源 的 亮度 分 布 图 像 可 以 通过 
测量 垂直 于 源 平面 上 的 点 对 之 间 电 场 的 互相 干 函 数 得 到 。 可 以 从 最 简单 的 没有 下 变频 、 没 有 延 时 补偿 
的 双 天 线 单元 接收 系统 进行 分 析 *] 。 

假设 s 为 天 线 指 向 射电 源 的 单位 方向 矢量 ， 两 天 线 的 基线 距离 为 了»， 其 中 方向 矢量 与 基线 法 线 方 
向 的 夹 角 为 90，V.,.(27f) 为 天 线 接收 系统 接收 的 来 自 方向 矢量 s 的 频率 为 的 单 频 信号 。 由 于 几何 距 
离 不 同 ， 射 电源 辐射 的 电磁 波 到 达 两 个 天 线 的 时 间 不 一 致 ， 产 生 的 时 间 差 为 


T, = 二 = jsingye， (1) 
C 


HP, 为 光速 ， 则 经 过 传输 链 路 到 达 相 关 融 之 前 的 输入 信号 可 以 表示 为 
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Vi (t) = Vicos [2af(t - 7,)], (2) 
V, (t) = V,cos 27fi, 
其 中 ，V, 、 妈 分 别 为 两 路 输入 信号 的 电压 幅度 ， 且 忽略 了 链 路 增益 。 相 关 器 由 一 个 电压 乘法 串 上 一 个 
时 间 平 均 (积分 ) 器 组 成 。 经 过 相关 器 后 滤 掉 高 频 成 分 分 量 ， 最 后 相关 器 的 输出 结果 为 
— [V,V,cos(2af7,)] LAL 
2 


cos (2mbsinb/àÀ ) , (3) 


其 中 ,入 为 观测 频率 对 应 的 波长 。 
当 射 电源 处 于 0=0 时 ,来 自 射电 源 的 信号 同时 到 达 接 收 天 线 ， 从 而 在 其 到 达 相 关 器 时 两 路 信和 号 
的 相位 差 为 0， 相 关 输 出 的 幅度 最 大 。 定 义 当 两 路 信号 到 达 相 关 融 时 相位 一 致 ( 即 输 出 幅度 最 大 时 ) 信 
号 的 来 波 方向 9 为 相位 中 心 ， 也 称 为 相位 参考 方向 。 在 实际 观测 中 ， 天 线 指向 偏离 目标 源 AO, MA 
是 处 于 相位 中 心 。 此 时 ， 考 虑 偏离 相位 中 心 Ab( Ab 很 小 ) 的 点 源 除 幅 度 外 的 相关 输出 结果 ， 表 达 式 
如 下 : 
cos(27fT,) = cos [2af [sin - A0) - "| = cos [2af (=) sinAdcos 0, |. (4) 


z= 在 以 上 计算 过 程 中 ， 忽 略 了 天 线 方向 办 的 影响 。 工 程 上 ， 天 线 接收 的 信号 强度 是 在 天 线 方向 图 范 
围 内 的 积分 。 假 设 辐射 源 的 亮度 分 布 为 2 ) ， 天 线 方向 图 为 4(9) ， 相 位 中 心 方向 为 (6 ) ， 天 线 指向 
为 (6o+Ab) ， 则 天 线 接收 到 的 信号 强度 为 


T V= AC, + A0) 1(0")d6", (5) 


其 中 ,，4(6,+Ab) 为 天 线 的 方向 图 ， 包 含 功 率 方 向 图 和 相位 方向 图 。 

a 2.2 测量 方法 

~N 干涉 仪 的 校准 和 测试 ， 通 常 采用 对 天 空中 的 点 源 进 行 观测 。 对 基于 相关 结果 测量 干涉 仪 天 线 相位 
> 方向 图 的 方法 ， 步 又 如 下 : 

N (1) 将 所 有 天 线 同时 指向 校准 源 ， 保 持 5 min 并 记录 此 过 程 中 的 可 见 度 函数 ; 

> (2) 保持 其 中 一 个 天 线 4 跟踪 校准 源 不 变 ， 将 其 余天 线 在 赤 经 方向 顺 时 针 逐 步 移 开 一 定 角 度 ， 保 
FF 1 min 并 记录 此 过 程 中 的 可 见 度 函数 ，; 

(3) 重 复 (1) 过 程 ， 然 后 保持 天 线 A 跟踪 校准 源 不 变 ， 将 其 余天 线 在 赤 经 方向 道 时 针 逐 步 移 开 一 
定 角度 ， 保持 1 min 并 记录 此 过 程 中 的 可 见 度 函 数 ，; 

(4) 重 复 (1) 过 程 ， 然 后 保持 天 线 A 跟踪 校准 源 不 变 ， 将 其 余天 线 在 赤 纬 方向 上 转 逐 步 移 开 一 定 
© 角度 ， 保 持 1 min 并 记录 此 过 程 中 的 可 见 度 函 数 ; 

(5) 重 复 (1) 过 程 ， 然 后 保持 天 线 4 跟踪 校准 源 不 变 ， 将 其 余天 线 在 赤 纬 方向 下 转 逐 步 移 开 一 定 
角度 ， 保 持 1 min 并 记录 此 过 程 中 的 可 见 度 函 数 ; 

(6) 根 据 可 见 度 函 数 计算 除 4 外 天 线 的 相位 变化 。 该 变化 即 是 由 于 天 线 指向 偏离 带 来 的 天 线 相位 
的 变化 ， 得 到 各 天 线 在 4 个 方向 的 相位 方向 图 。 

对 于 MUSER-I， 指 向 精度 和 跟踪 精度 9'。 位 于 赤道 上 方 的 同步 轨道 卫星 是 非常 合适 的 校准 源 ， 
根据 MUSER-I 的 观测 频率 ， 风 云 2G 卫星 是 一 个 合适 的 校准 源 。 风 云 2 G 卫星 是 我 国 的 一 颗 地 球 静 止 
气象 卫星 ， 位 于 东经 104. 9 赤道 的 上 空 ， 其 信 标 频率 为 1.702 5 GHz。 由 于 卫星 信号 比较 强 ， 并 且 在 
该 频 点 上 干扰 信号 非常 小 ， 所 以 在 测试 过 程 中 可 以 忽略 干扰 。 

2015 年 10 月 28 日 20 点 ,按照 上 述 步 又 利用 风云 2G 卫星 (时 角 : 12. 10” ， 赤 纬 : -4. 86° ) 测 量 
了 MUSER-I 天 线 的 相位 方向 图 。 

2.3 测量 结果 与 分 析 

通过 2. 2 中 记录 的 可 见 度 函数 可 以 计算 相位 值 ， 如 图 2 和 图 3。 其 中 图 2 表示 天 线 指向 在 赤 经 方 
向 偏离 相位 中 心 1° 时 的 相位 差 ， 图 3 表示 天 线 指向 在 赤 纬 方向 偏离 相位 中 心 +1° 时 的 相位 差 。 由 于 
太阳 的 尺度 大 概 为 30"， 所 以 重点 分 析 天 线 在 相位 中 心 上 zl 范围 内 的 相位 方向 图 。 根 据 图 2 和 图 3 在 
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+1 x+1 的 范围 内 ， 天 线 的 相位 与 相位 中 心 的 相位 差 基本 在 5 范围 内 ， 而 在 太阳 的 成 像 范围 +30' x + 
30' 内 ， 天 线 的 相位 与 相位 中 心 的 相位 差 基 本 在 2 范围 内 。 

在 测量 过 程 中 由 于 卫星 的 运动 会 引起 日 像 仪 输出 相位 的 变化 ， 在 数据 处 理 过 程 中 根据 源 的 位 置 变 
化 和 UV 变化 的 关系 ， 算 出 由 于 位 置 变化 引起 的 相位 变化 ， 从 测量 结果 中 扣除 。 

由 于 太阳 的 展 源 特性 ， 成 图 主要 取决 于 天 线 主 匆 的 幅度 方向 图 和 相位 方向 图 的 平坦 情况 ， 而 天 绑 
的 偏离 相对 于 幅度 误差 会 引起 大 的 相位 误差 。 因 此 在 分 析 过 程 中 可 以 忽略 幅度 误差 而 重点 分 析 相 位 误 
差 对 成 图 质量 的 影响 。 

相位 误差 对 图 像 的 影响 取决 于 其 对 可 见 度 函数 的 影响 。 通 常 基于 天 线 的 误差 和 基于 基线 的 误差 是 
可 见 度 函数 测量 中 最 普遍 的 两 种 误差 ,天 线 相 位 方向 图 便 属 于 基于 天 线 的 误差 。 
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图 2 天 线 指向 在 赤 经 方向 偏离 相位 中 心 上 1" 时 的 相位 差 


Fig. 2 Measured phase difference when antenna pointing deviates + 1° in the right ascension direction 
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图 3 天 线 指向 在 杰 纬 方向 偏离 相位 中 心 +1 时 的 相位 差 


Fig.3 Measured phase difference when antenna pointing deviates + 1° in the declination direction 


对 于 明 安 图 射电 日 像 仪 考虑 简单 的 位 于 相位 中 心 单位 流量 密度 的 点 源 ， 使 用 NN 个 天 线 的 干涉 仪 
观测 ， 在 快照 模式 下 可 见 度 函数 为 


Vu) = (1 +e)6(u - uy) e”, (6) 
其 中 ，s 为 可 见 度 函 数 中 的 幅度 误差 ， 为 可 见 度 函数 中 的 相位 误差 ;wu, 为 干涉 仪 的 基线 。 
干涉 仪 的 输出 图 像 可 以 表示 为 
I0) = [VGu) ex™qu. (7) 


因为 干涉 仪 的 输出 是 WCN-1)7Z2 个 离散 的 分 量 ， 所 以 积分 可 以 变 为 求 和 计算 。 对 于 每 个 不 存在 误差 
的 基线 ， 对 图 像 的 作用 为 2cos(2mz1) 。 但 是 对 于 存在 误差 的 基线 ， 这 一 分 量变 为 2cos(2mrxl- 中 ) 二 2 
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[cos(27u 1) +qsin(2mwuol)]。 因 此 ,干涉 仪 最 终 输 出 的 图 像 可 以 表示 为 


N(W=1)/2 
I(l) = 2dbsin(2Tu,l) + 2 > cos(27u, 1) , (8) 
此 时 干涉 阵 的 点 扩散 函数 可 以 表示 为 
N(N-1)/2 
B(l)=2 > cos(27u,l). (9) 


通过 对 输出 图 像 做 洁 化 处 理 即 去 卷 积 ， 可 以 得 到 误差 图 像 
R) = I(LI) - BO) = 2dsin(271u, 0) , (10) 
这 一 误差 图 像 是 一 个 正弦 函数 ， 其 幅度 为 26$， 有 效 值 为 V24$。 
图 像 的 动态 范围 可 以 定义 为 图 像 中 辐射 强度 最 大 值 与 一 个 无 源 小 区 域内 的 均 方 根 的 比值 。 通 过 上 
述 公 式 的 推导 可 以 得 到 : 当 干 涉 仪 由 NN 个 天 线 组 成 ， 对 应 有 N(N-1)/2 个 基线 ， 可 得 到 N(N-1)/2 
个 可 见 度 函 数 ， 假 定 其 中 的 一 个 基线 具有 相位 误差 $9， 此 时 干涉 仪 输出 的 动态 范围 为 


Peak N(N-1) M 
D = 一 = = (N>1Ħf). (11) 
= Noise 26 26 
A 现实 中 误差 很 少 只 存在 于 一 个 基线 上 ， 所 以 假定 所 有 的 基线 上 存在 一 个 随机 的 相位 误差 ， 并 且 每 个 基 
is 线 上 的 误差 是 独立 的 ， 则 图 像 的 动态 范围 随 着 比例 因子 VNCN - 1) 72 下 降 ， 即 
a N(N - 1) J/N(N-1) N 
© D= = =—. (12) 
ia J26/NC(N - 1)72 由 由 


如 果 相 位 误差 是 基于 天 线 的， 误差 影响 与 该 天 线 相关 的 所 有 基线 ， 所 以 对 图 像 造成 N-1 个 误差 。 假 
设 这 N-1 个 误差 是 随机 的 ， 则 图 像 的 动态 范围 可 以 表示 为 
N _ N(N-1) _NJSCN=T) ~ 应 N 
> BeN- 6 Tg 
同样 现实 中 所 有 天 线 有 基于 天 线 的 相位 误差 ， 此 时 图 像 的 动态 范围 随 着 比例 因子 WW 下 降 ， 可 以 表示 为 


re ea: (14) 
$ J2b 


(13) 


D 


© 因此 ， 对 于 日 像 仪 低频 阵 ， 天 线 数目 为 40 面 ， 技 术 指标 要 求 动态 范围 25 dB(320 : 1) ， 根 据 公式 计 
算 日 像 仪 的 容许 相位 为 5 。 可 见 在 太阳 成 像 范围 内 日 像 仪 天 线 的 相位 差 在 2* 范 围 内 满足 要 求 。 


3 结 论 


文中 根据 日 像 仪 的 馈 源 设计 和 综合 孔径 原理 ， 针 对 明 安 图 射电 频谱 日 像 仪 天 线 数目 多 ， 且 为 户外 

环境 ， 根 据 天 文 观测 须 经 常 测试 天 线性 能 的 特点 ， 给 出 了 基于 相关 结果 测量 日 像 仪 天 线 相 位 方向 图 的 

方法 ， 由 图 2 和 图 3 可 知 明 安 图 射电 频谱 日 像 仪 天 线 的 相位 差 满 足 太 阳 成 像 动态 范围 25 dB 的 技术 
旨 标 要 求 。 
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Abstract: The MUSER ( Mingantu Ultrawide SpEctral Radioheliograph ) is a new generation of solar 


dedicated radio telescope for imaging the Sun at high time, high spatial and high frequency resolutions 


simultaneously. It is an aperture-synthesis telescope, so the amplitude and phase are the key factors in the 


final imaging. The phase pattern is an important factor affecting the array’ s amplitude and phase. In this 


paper, allowing for the techniques of MUSER, the method by measuring correlation is adopted to measure 


phase pattern of feeds for MUSER based on aperture synthesis principle. Since MUSER has numbers of 


antennas and located in harsh environments and according to the demand of astronomical observations, the 


performance of the antennas should be measured frequently. This method simplifies the process of phase pattern 


measurement of the interferometric array and it can obtain accurate results efficiently. The measuring process 


and analysis of MUSER-I are presented in detail in this paper. In addition, to support future scientific 


observation, the phase pattern’ s effect on the image is further studied. 


Key words; Phase pattern; Correlation; Mingantu Ultrawide SpEctral Radioheliograph 


